Spontaneous emission interference in negative-refractive-index waveguides 
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The spontaneous decay of a V-lype three-level atom placed in a negative-refractive-index waveguide is ana- 
lyzed. We find that in thin waveguides, highly efficient surface guided modes are supported, which do not occur 
in positive index waveguides. In addition, at low absorption, the mode density and thus spontaneous emission 
into particular regular guided modes is enhanced by several orders of magnitude as compared to regular dielec- 
tric waveguides. The asymmetries between emission into the different modes and the enhancement of particular 
guided modes allow to induce strong spontaneous-emission interference between transitions with orthogonal 
transition dipole moments. 

PACS numbers: 78.20.Bh,78.20.Ci,42.50.Gy,42.50.Pq 



I. INTRODUCTION 

Spontaneous emission of atoms is not an immutable prop- 
erty, but can be altered essentially via two different mecha- 
nisms. One mechanism involves modification of the internal 
dynamics. For example, quantum interference among dif- 
ferent decay channels of the atom can be established, such 
that spontaneous emission is modified to a great extend'-^. A 
particular class of quantum interference schemes that has re- 
ceived much theoretical attention in the literature is based on 
so-called spontaneously generated coherences (SGC). These 
coherences have for example been shown to lead to quench- 
ing of spontaneous emission^, narrow spectral lines'*, phase- 
dependent line shapes^, and rapid phase control of collec- 
tive population dynamics^. Despite the large theoretical in- 
terest, there is no experimental proof of this type of SGC in 
atomic systems due to the lack of appropriate candidate sys- 
tems. The reason for this is that SGC based interference re- 
quires the presence of near-degenerate atomic transitions with 
near (anti-)parallel dipole moments sharing a common atomic 
state, which does not occur in real atoms^. To circumvent this 
in atoms, schemes to simulate interference or to induce inter- 
ference by external driving fields have been considered'-^, but 
so far without clear experimental implementation. An exper- 
imental observation of SGC in molecules** could not be con- 
firmed in a repetition of the experiment^. SGC have been ob- 
served, however, in artificial quantum systems, where a suit- 
able level scheme can be designed'". Also, it has been recog- 
nized that a more general form of SGC occurs between near- 
degenerate atomic transitions with near (anti-)parallel dipole 
moments that do not share a common state. This generalized 
SGC may lead to measurable effects in realistic atomic sys- 
tems". These generalized SGC are also a common interpre- 
tation for the occurrence of electromagnetically induced ab- 
sorptioni^. 

The second major mechanism for a modification of spon- 
taneous emission is to modify the electromagnetic boundary 
conditions surrounding the atom, such as in cavities^^. Re- 
cently, is has been shown that a favorable modification of the 
boundary conditions is possible in media having a negative 
refractive index (NRI). NRI currently receive a lot of atten- 
tion because of both experimental demonstrations and poten- 



tial applications'^. For example, NRI material allows to real- 
ize superlenses which, in principle, can achieve arbitrary sub- 
wavelength resolution'^. Experimentally, the NRI materials 
have been realized over a wide range of frequencies, from the 
microwave up to the optical range""''^. It has been shown '-^ 
that a single layer of a NRI material has broadband omnidi- 
rectional reflection properties. In the range between the elec- 
tric plasma frequency and the magnetic plasma frequency, the 
refractive index is close to zero and the NRI materials reflect 
radiation for angles of incidence and polarization with reflec- 
tivity of ^ 0.99. By exploiting the reflective phase proper- 
ties of the NRI materials, an all-dimensional subwavelength 
resonator" has been designed and fabricated. In the field of 
quantum optics, two atoms placed at the foci of a perfect lens 
formed by a NRI slab exhibit perfect sub- and supeiTadiance 
over macroscopic distances^". Such focussing and phase com- 
pensation can also be used to induce quantum interference^!. 
In a NRI slab waveguide^^, surface guided modes may exist 
for imaginary transverse wave numbers, with power concen- 
trated at the interfaces rather than inside the slab. In an air 
waveguide with NRI cladding-^, both transverse electric and 
magnetic modes can be supported with low losses. The exis- 
tence of the transverse electric (TE) surface guided modes in 
the NRI slab waveguide, where the permeability changes its 
sign at the interface, resemble the transverse magnetic (TM) 
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FIG. 1: (Color online) Schematic setup for spontaneous decay in a 
waveguide with negative index of refraction material, (a) The waveg- 
uide consists of an inner layer with material parameters es, fi^, and 
thickness da. It is surrounded by infinite outer layers with material 
parameters e^, fii (i £ {1, 2}). In the middle layer, an atom is em- 
bedded at position ro. (b) Electronic structure of the embedded atom, 
which is a three-level atom in V-type configuration. 
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surface-plasmon-polariton modes at the interface between the 
metal and the dielectric, where the permittivity changes its 
sign^'*. Plasmonic -based nanophotonic devices have attracted 
much interest from the quantum optics community for their 
use in quantum information processings^. A technique which 
enables strong, coherent coupling between individual optical 
emitters and guided plasmon excitations in conducting nanos- 
tructure at optical frequencies has been proposed^''. In view 
of these remarkable properties of NRI structures, the question 
arises, how spontaneous emission is modified in such a sur- 
rounding. 

Therefore, here we study the spontaneous emission (SE) of 
an atom embedded in the middle layer of NRI slab waveg- 
uides. The atom is modelled as a V-type three-level atom with 
orthogonal dipole moments on the dipole-allowed transitions 
as it is the case in real atoms. The NRI layer is described using 
a Drude-Lorentz model including dispersion and absorption. 
The different contributions to SE are first classified and inter- 
preted using approximate analytical results. Then, we verify 
our results via numerical calculations. We find that at low ab- 
sorption, the SE rate into particular waveguide modes can be 
several orders of magnitude larger than the free space rate due 
to a strong enhancement of the mode density. Further, in thin 
waveguides, a dominant contribution to SE arises from surface 
guided modes, which do not occur in regular positive index 
waveguides. At high absorption, the strong enhancement of 
waveguide modes is reduced, while the surface guided mode 
contribution remains several orders of magnitude higher than 
free space decay. As our main result, we show that the asym- 
metric mode structure and the strong enhancement of partic- 
ular modes can be used to engineer spontaneous emission in- 
terference between the two transitions in the embedded atom. 
Thus we find that the NRI waveguide allows to induce and 
control near-perfect quantum interference in realistic atomic 
level schemes. 



dipole strength, chosen to be real. We denote the dielec- 
tric permittivity and the magnetic permeability of the layer 
cladding the atom as e3{uj) and ij.3{uj), which we obtain from 
a Drude-Lorentz model: 
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Here, ajp^ (ujpm), ^Te i^Tm) and 7^ (7„i) are electric (mag- 
netic) coupling constant, medium oscillation frequency and 
linewidth. The permittivity and permeability of the upper 
[lower] layer, which extends to infinity in the positive [neg- 
ative] z direction, are e2(t^) and /i2(w) [61(0;) and ^i(w)]. 
For these, we assume low absorption and choose ej, fij (j G 
{1, 2}) as real22i22i21. ^3 is the thickness of the middle layer. 
The interaction picture Master equation for the density matrix 
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Here, r„ are spontaneous emission rate 
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with polarization components and Tnz- The terms involv- 
ing Ki and K2 in Eq. (|4]i are responsible for quantum interfer- 
ence between the two SE channels |1) |3) and |2) |3), 
with K„ given by 



n. THE MODEL 
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We consider a V-type three-level atom embedded in 
the middle layer of a three-layer waveguide as shown in 
Fig. [1] The upper levels could be Zeeman sublevels 
|1) = |j = l, 771=1) and |2) = |j = l,m = — 1) with energies 
huji, huo2, and we set the ground state |3) = |j = 0,777 = 0) 
energy to zero in the following. The y direction is the quanti- 
zation axis (e.g., by applying a weak static magnetic), and the 
direction normal to the layer interfaces the z axis. Then, the 
atomic dipole moment operator is given by 



where 



d = d(Ai3ei + ^2362) + H.c. 
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and A. 



I*) 01 (*>i £ {lj2,3}) are the atomic transition 
(i 7^ i) and population (i — j) operators, e^. (fc e {x, y, z}) 
are the normalized Cartesian basis vectors, and d is the atomic 



— 1 < Ki,K2 < 1 describes the degree of interference, and 
To = (a^O: yo, zq)'^ is the position of the atom in the middle 
layer Gim = Im [G(ro, rp, ujn)] is the imaginary part of the 
electromagnetic Green tensor G(ro, Tq, lu) given by"** 
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Here, 



Xi"^ = (1 ± r|ie2*^^^°)(l ± rl2e^'^'^'''-''>^) 

l2 2/2 

r]j = £j(w)^j(w), 



(7) 
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with i G {1, 2, 3}. The parameter k is the magnitude of the 
vector k = [kx, kyf^ , the conserved component of the wave 
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FIG. 2: (Color online) (a) Radiation modes, with total rate Tr (red 
solid line) and contribution T^^ of the p-mode for z-polarization 
(blue dashed), (b) Waveguide modes, with total rate Vg (red solid 
line) and contribution Fg^ of the p-mode for z-polarization (blue 
dashed). The atom is located at z'^ = 0.25d3, and the material pa- 



rameters are ei = £2 = Mi = ^2 = 1-0, Upe = i^^pm ~ 1.32[i;o, 

LOTe = LOTm = LOQ, (^a = 1.09Wo, and Je = 'fm — W'^'^LdQ- 

Here, ujo is a scaling parameter. Then, ^^{uja) ~ e3{uja) ~ 
-1.99 + 1.73 -IQ-^i. 



vector, which is parallel to the interfaces of the layers. The 
Pj [j € {1, 2,3}) are the magnitude of the z component of 
the wave vector in the j — th layer, whose definition depends 
on the refraction index of the j — th layer and the value of 
h^. If Re(A;) < Re(fcj), that is, the corresponding wave in 
the j — th layer is a propagating one, then [3j is expressed as 
(3j — {k^ — k^y/^ if the j — th layer has positive refraction 

index, and Pj — — (fcj — k^y^^ when the j — th if the layer 
is a left-handed material. Here, Re(x) is the real part of x. 

On the other hand, if Re(fc) > Re{kj) corresponding to an 
evanescent wave in the j — th layer, then fdj — i{k^ — kjY^^ 
independent of the type of the j ~ th layer Following Ref.^^, 
we denote the electric field of the TM [TE] wave by the index 
p [s]. The functions Dg^ (q e {p, s}) are defined as 



— '3l'32'^ 



where r|]^ and are reflection coefficients given by 



(9) 
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From Eq. Q we can see that the inhomogeneity of the 
medium along the z axis leads to a spatially asymmetric Green 
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FIG. 3: (Color online) Decay rates into surface guided modes, with 
total decay rate contribution Fs and contribution Ff ^ of the p-mode 
for z-polarization. The parameters are as in Fig.|2] 



tensor, as the z-component is different from those in the x — y 
plane. For notational simplicity, in the following and k 
will be re-scaled by ui/c, i.e, ^^/{llj/c) and k/{u!/c) will 
be replaced by P3 and k. As usual, the permittivity and 
permeability coefficients of the three layers are assumed to 
obey Re (773) > Re (771) > Re (772)- Then the electromag- 
netic modes in this structure can be classified into radiation 
modes with < Re{k^) < Re{ri2), substrate modes with 



Re{m) 



i?e(fc2) 
i?e(fc2) 



< i?e(?7i), regular guided modes with 

< Re{r]3), and evanescent modes with 
Re{k^) > Re{ri3). As quantum interference requires near- 
degenerate transition frequencies, we also assume uji uj2 — 
Lo, SO that Fi « r2 = r and ni k, H2 = n. 



III. RESULTS 

We assume that both the upper and the lower layers are di- 
electric media with positive refraction index, and the middle 
layer surrounding the atom is left-handed material. This struc- 
ture is a so-called negative-refraction-index (NRI) waveg- 
uide-^. The contributions of the radiation and the substrate 
modes to the SE rate can be obtained by integrating Eq. (|7]) 
numerically. An example is shown in Fig. |2|a). For the fol- 
lowing analytical considerations, we assume £3 and /i3 to be 
constant and real, i.e., we neglect absorption and dispersion. 
But our numerical results include both absorption and disper- 
sion. 



A. Waveguide modes 

The regular guided modes have complex reflection coef- 
ficients with modulus 1. In this region, the imaginary parts 
of the integrands in Eq. ^ are zero, apart from resonances 
when Z)q3 = (q G {p, s}) is satisfied. For example, the SE 
rate of the z— component of the atomic dipole moment into 
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p— polarized guided modes is 
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in which Zg = zquj/c, d'^ — d^uj/c, and 70 represents the SE 
rate of the atom in free space. Further, 



rfa = exp(-2i0P2) 



withO < (t>l2 < 7r/2, and 



Xi 



£3A<3 - £^^J'^ 



(12) 



(13) 



The parameters kg^^ (m e {1,2,..., m-^J^^^}) represent the 
wave numbers of the mth p— polarized guided modes in the 
X — y plane, which are the real roots of equation Dp^ = 
within the region of y/rfi < k < ^/fj^. The number of modes 
is mf^^^, which depends on the thickness of the middle layer, 
the mode polarization, and the material parameters. Evidently, 
each guided mode corresponds to a standing wave in this 
structure, and the SE rate depends on the position of the atom. 
With increasing thickness ^3, the contribution of the guided 
modes to the SE rate exhibits a sharp cusp at the appearance 
of an extra mode. For the symmetric case ei = £2, /ii — /i2 
and zq = ^3/2, there is no contribution to the decay rate from 
the p— polarized modes whose nodes are exactly coincident 
with the atomic location, such that (S^ds + "^4^2 = (^^^^ + l)7r 
{n E {0, 1,2,... }). Only modes with an anti-node coincid- 
ing with the atomic position contribute to the decay rate, since 
then Psds + 2^32 — 2mT. 

Full numerical results for the guided modes are shown in 
Fig.|2b). Note that guided modes may already occur for very 
thin, even sub-wavelength, layers, as the Goos-Hanchen phase 
shifts (GHPS) at the interfaces between negative and positive 
layers enhance the phase change induced by the optical path^^ . 
The second and the third terms in the denominator of Eq. (fTTl i 
arising from the GHPS are negative, while the first term due to 
the optical path is positive. Via this cancelling, in z-direction 
the structure effectively acts as a planar cavity with tiny length 
and large mode density. In Fig.|2l the decay rate is enhanced 
by about 4 orders of magnitude compared to the free space 
decay rate for specific waveguide geometries. This is in sharp 
contrast to the regular dielectric waveguide with £3 , ^3 > 0, 
where all three terms in the denominator are positive, such that 
the guided modes have only a small density of modes. Thus 
in a suitable NRI waveguide, the atom spontaneously emits 
photons into few guided modes with large amplitudes. For an 
atom placed close to a nano-structure, spontaneous emission 
can be greatly enhanced due to couplings to electronic quasi- 
particle surface excitations^i. But here, the enhancement of 
some guided modes arises from a variation of the summation 
of the optical path and the phase changes at the interfaces be- 
tween negative and positive refraction index materials. Thus 
these interfaces lead to the strong enhancement. 
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FIG. 4: (Color online) The influence of dispersion and absorption on 
the results. The red solid line shows the frequency dependence of 
the total decay rate for parameters as in Fig. (2] but with dg set to the 
first peak in . The blue dashed line is plotted for same parameters 
but includes small absorption 7 = 7e = 7m ~ W^^loo. The black 
horizontal dotted lines indicate the peak amplitudes of about 7IOO70 
and 7IO70, respectively, of the two curves, verifying the scaling of 
the peak amplitudes with 1 / ^7. 



B. Surface guided modes 

In regular dielectric waveguides, evanescent modes do not 
contribute to the SE rate, as they cannot propagate out. But if 
one or two of the layers are made from NRI, surface guided 
modes with fc^ > 7^3 can exist, whose wave vectors obey the 
resonance condition Dq^ = in the Green's function Eq. 
Different from the regular positive index dielectric waveguide, 
where the magnitudes of the reflection coefficients and 
r|]^ are smaller than 1, in the NRI waveguide, the modulus of 



and r|j^ can be larger than 1. Thus real roots fci™-* (m g 
{1,2,...}) can be found, corresponding to the wave vectors 
of the surface guided modes. The contribution of these guided 
modes to the z— component of the dipole moment is 
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(14) 



k=k), 



Here we have set = ~ exp(2032) with (j>^2 ^ 0- 
In Eq. (fT4b . the mode functions are hyperbolic instead of 
the standing wave functions for the regular guided modes 
Eq.dnj. 

By numerical and analytical inspection, we find that there 
are two kinds of surface guided modes. The first has a wave 
vector k^ close to 773, while the other one has k^ ^ rj^. It can 
be proven that the surface guided modes with fc^ ^ 773 can 
only exist when the thickness ^3 of the middle layer is very 
thin. As an example, we estimate the condition for this kind 
of surface guided modes for the symmetric NRI waveguide 
(ei = £2, /^i = A'2). For the p— polarized surface guided 
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FIG. 5: (Color online) Spontaneous emission into strongly absorbing 
NRI waveguides, (a) Radiation modes, with total rate Vr (red solid 
line) and contribution T^^ of the p-mode for z-polarization (blue 
dashed), (b) Waveguide modes, with total rate Tg (red solid line) 
and contribution F^^ of the p-mode for ^-polarization (blue dashed). 
The atom is located at Zq = 0.5^3, and the material parameters are 

£l = £2 = A*! — M2 = 1.0, LUpe ~ 1.25tJo, UJpm — l.lSOoJo, 

lotc ~ i^Tm = <^o, i^a = l.OStJo, and 7e — 7m ~ W^'^ujo. Here, 
LOO is a scaling parameter. Then, e^luja) ~ —2.38 + 2.19 ■ lO^^i 



and ^i3{(^a) 



-1.48 + 1.61 • 10' 



modes, the condition can be approximated as 
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Therefore polarized surface modes with fc^ ^ £3/^3 may 
exist only if the thickness obeys ^ In | (ei — £3)/(£i + 
£3)|/'73, i-C-, 0^3 is much smaller than one wavelength. For 
these surface modes, Eq. (fT4b reduces to 
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We can see that if Zq ^ 0^3/2, then the decay rate can be large, 
as shown in Fig. [3] as this mode only survives for thin ^3 with 
a large amplitude. This large enhancement of surface guided 
mode excitation enables dipole emission to be preferentially 
coupled to the surface guided modes, which may be applied 
in creating well-guided light sources at the nanoscale. 
Regarding the second type of surface guided modes with 
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FIG. 6: (Color online) Decay rates into surface guided modes. The 
total decay rate Fs has contributions only from Ff^. of the p-mode 
for a;-polarization. The parameters are as in Fig.jS] 



close to ?73, it can be found that these modes may exist if 
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This is the reason for the sharp cutoff visible for the solid 
curves in Fig. [3] The decay rate induced by these p— polarized 
surface guided modes can be approximated as 
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In thin waveguides, these contributions are small as compared 
to those of Eq. ( fT6b . due to the difference of the mode ampli- 
tudes. 

Thus for thin NRI waveguides, the atom can emit strong 
evanescent fields with large wave number fc^ ^ VslJ-s in 
the X — y plane. This result is similar to experimental evi- 
dence employing a pristine silver film with natural roughness 
as a NRI slab-'^, where the transmission of evanescent waves 
rapidly grows with the NRI film thickness up to a thickness of 
about A/ 10. 



C. Dispersion and absorption 

With absorption and dispersion, the above equations are 
more complex. Then, the sharp peaks due to the guided modes 
become Lorentzians, with widths depending on the imaginary 
part of the refraction index. For small 7 = 7e = 7m, their 
amplitudes are proportional to This scaling generalizes 

our results to moderate absorption strengths. An example is 
shown in Fig.|4] This figure shows a magnification on the first 
peak at dg w 3 in the waveguide decay F^^ shown in Fig.|2l It 
can be seen that with increasing absorption, the naiTow peak 
is broadened such that the peak amplitude is reduced, in ac- 
cordance with our scaling law. In contrast to Fig.|2l in Fig.|4] 
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the frequency-dependence of the resuhs is shown in order to 
analyze the dispersive properties of the NRI waveguide. 

We now turn to the case of high absorption. An example for 
decay into radiation modes in a strongly absorbing waveguide 
is shown in Fig. HJa). The corresponding results for waveg- 
uide modes and special guided modes are shown in Fig.|3b) 
and Fig. |6l It can be seen that while the waveguide modes are 
strongly reduced to peak decay rates of order 70, the radia- 
tion modes and special guided modes are much less affected 
by the absorption. Speaking pictorially, the absorption leads 
to a washing out of narrow structures found in the waveguide 
modes, such that their amplitude is strongly reduced with in- 
creasing absorption. In contrast, the structures in the radiation 
and special guided modes are akeady rather broad at low ab- 
sorption such that they do not change much towards higher 
absorption. 



D. Quantum interference 

Finally, we evaluate the cross-coupling terms k„ in Eq. (|4|, 
which are responsible for SE interference. The strong en- 
hancement of decay into particular guided modes at low ab- 
sorption can be used to generate controllable, near-perfect 
quantum interference, i.e. k sa ±1, as shown in Fig. |3a). 
If both transitions mainly interact with the same strongly en- 
hanced guided mode, then they couple even though they have 
orthogonal transition dipole moments. For example, if ^3 is 
chosen close to the birth of a s— polarized mode, but away 
from p-polarized mode maxima, then the decay rate is dom- 
inated by induced by the s— polarized modes to the x- 
component of the atomic dipole. It follows from Eq. ^ that 
then K sa — 1, i.e., strong quantum interference similar to 
that of a hypothetical atom having two near-degenerate tran- 
sitions with near-antiparallel dipole matrix elements in free 
space '^ii^. A second example is for atoms located at the 
center of the NRI layer, with set at the birth of a p-polarized 



(a) 



mode. If /Jgdg + 20^ — {2n + l)7r, then the decay rate 
is dominated by F^^ from the component of the atomic 
dipole, as the atom is at the node of the z-modes. But for 
/Sadjj + 2(/)^ ~ 2mr, it is dominated by Fg^, as the atom is at 
the node of the x-modes. These situations yield strong quan- 
tum interference with k « ±1, similar to near (anti-) paral- 
lel dipole moments in free space. Both cases are shown in 
Fig. |2f a). The extremal values of k occur at peaks in the SE 
rate due to the enhanced mode density, and thus d'^ controls 
K between —1 and 1. In Fig.|2la), we have used 63 ^ /Z3 to 
enable the first mechanisms for extremal k, and Zq = ^3/2 to 
allow for the second mechanism. Thus NRI waveguide struc- 
tures allow to effectively induce spontaneous-emission quan- 
tum interference with realistic atomic level structures. 

Interestingly, spontaneous emission interference can al- 
ready be obtained at high absorption, even though in this 
case the narrow peaks of large amplitude are absent from the 
waveguide mode spectrum. An example is shown in Fig.|7|b) 
for the parameters of Fig.|5] In particular at small thicknesses 
1^3 of the waveguide, a large degree of interference is achieved 
(k Ri —1). The reason for this is the dominant contribution of 
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FIG. 7: (Color online) Quantum interference induced by a negative- 
refraction-index waveguide structure, characterized by the interfer- 
ence strength n. (a) Parameters as in Fig. [2] except for z'q = 



-1.20 + 1.27 ■ lO^-'i]. (b) Strong ab- 
i -2.38 + 



sorption case. Parameters are as in Fig.jS] with £3 (a 
2.19 ■ lO^^z and fisicoa) ^ -1.48 + 1.61 • 10" ' 



the special guided modes. For the particular setup chosen in 
this figure, only special waveguides in the p — x mode are 
excited. But also at larger thicknesses, interference of alter- 
nating sign in k is achieved, even though with |k| smaller than 
unity. This has to be compared, however, to a value of k = 
for atoms in free space. 



rV. SUMMARY 

In conclusion, we have investigated the spontaneous decay 
of a three-level V-type atom placed in a three-layer waveguide 
with negative-refraction-index material as middle layer We 
have found that spontaneous emission into particular guided 
modes can be greatly enhanced in materials with low losses. 
Due to a large mode density, spontaneous decay rates into 
these modes can be increased by several orders of magni- 
tude as compared to regular dielectric waveguides. Both at 
low and high absorption, NRI waveguides support additional 
surface guided modes, which do not occur in positive index 
waveguides. The modes also offer a strong enhancement of 
spontaneous decay compared to the free space case, and are 
especially effective in thin waveguides. The specific proper- 
ties of the waveguide as well as the position of the atom in the 
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waveguide enable one to modify the emission into the differ- 
ent modes to a great extend. We have shown that this control 
can be used to achieve asymmetries in the decay into modes 
with different polarizations. This feature allows to induce 
strong spontaneous emission interference in realistic atomic 
level schemes, which do not exhibit such interference in free 
space. 

In addition to NRI waveguides, we also analyzed air-like 
waveguides with NRI cladding^-'. In these structures, the mid- 
dle layer is air-like, while the upper and the lower layer are 
made of NRI material. We obtained similar results as for the 
NRI waveguides. In particular, also in this case the density of 
the guided modes and thus SE can be very large if the absorp- 
tion is low. 
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